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Abstract 

The probability distribution for the relative return of a portfolio constructed from 
a subset n of the assets from a benchmark, consisting of N assets whose returns are 
multivariate normal, is completely characterized by its tracking error. However, if the 
benchmark asset returns are not multivariate normal then higher moments of the proba- 
bility distribution for the portfolio's relative return are not related to its tracking error. 
We discuss the convergence of generalized tracking error measures as the size of the 
subset of benchmark assets increases. Assuming that the joint probability distribution 
for the returns of the assets is symmetric under their permutations we show that in- 
creasing n makes these generalized tracking errors small (even though n << N). For 
n >> 1 the probability distribution for the portfolio's relative return is approximately 
symmetric and strongly peaked about the origin. The results of this paper generalize 
the conclusions of Dynkin et. al. (2002) to more general underlying asset distributions. 



The excess return of portfolios managed against a benchmark can be attributed to active 
management in the selection of different weights from the benchmark, or from the selection 
of out-of-benchmark securities. For a portfolio that attempts to replicate a benchmark's 
risk and performance with a subset of the benchmark securities under a portfolio repli- 
cation approach, it is usually assumed that minimizing the tracking error of the portfolio 
with respect to the benchmark is equivalent to minimizing the underperformance or out- 
performance of the portfolio's relative return. In principle, this replication approach suffers 
from the fundamental problem that the tracking error as a measure of risk relative to the 
benchmark fails in its essential purpose if the underlying asset returns are not joint normal, 
i.e. when it is needed the most. In other words, if the asset returns are correlated and far 
from multivariate normal then it is possible that a tracking error optimized portfolio will 
not follow the benchmark as well as reasoning, based on the assumption that the portfo- 
lio's relative return is normal, would lead one to believe. But returns on corporate bonds 
subject to default risk or downgrade risk are not multivariate normal even if the underlying 
corporate assets are assumed to be since default of a company only occurs if the value 
of its corporate assets fluctuate below its liabilities^]. Furthermore, one does not have the 
luxury of increasing the number of assets in the portfolio ad infinitum since the addition 
of each new asset increases transactions costs and makes the portfolio less liquid. Thus it 
is important to consider, for a wide range of asset probability distributions, how the total 
risk of a portfolio relative to the benchmark saturates as more assets from the benchmark 
are kept. 

Consider a benchmark portfolio with N assets that return fj, i = 1, ... TV. Assuming 
the value of the benchmark portfolio is equally divided amongst these assets the random 
variable for the benchmark portfolio return is 
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Suppose a portfolio V consists of equal amounts of n of these assets which, without loss 
of generality we choose to be i = 1, . . . , n. Then the random variable for the return of V 
relative to the benchmark is 
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where, 
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when i = 1, . . . , n and 

when i = n + 1, . . . , N. The joint probability distribution for the asset returns can be 
characterized by the expected values, 

4™] im = E[r il ...hJ. (5) 



1 This is usually described using a first passage time model. See, Black and Cox (1976), Longstaff and 
Schwartz (1995), Leland and Toft (1996), etc. 
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We assume the joint probability distribution for the asset returns is symmetric under 
their permutations. This should make it as easy as possible for V to track the benchmark 
while still allowing for fluctuations in 7"s relative return. More realistic situations where 
the returns behave differently will not track the benchmark as well. With this assumption 
the expected asset returns are all equal, r- 1 ^ = \ the diagonal elements of the asset 



return covariance matrix are all the same, r\ 
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and the off-diagonal elements of the 



asset return covariance matrixF] are also independent of which asset is being considered, 
for i ^ j. Note however we do not demand that the off-diagonal elements be 
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equal to the diagonal ones since this would correspond to a correlation matrix that has 
one in all its elements. The generalization of this to higher values of m is straightforward. 
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-^2 f° r * 7^ 3 an d f^tjl = r u3 for i ^ j ^ /c, etc. 



Explicitly, r\f] 

In the presence of correlations the central limit theorem cannot be used to argue that 
for N >> 1 the benchmark return R becomes normally distributed or that for n >> 1 
Vs return becomes normally distributed. It is sometimes argued that for large portfolios 
correlations should average out leaving one in a situation where effectively the central limit 
theorem is applicable. However, we consider that to be very unlikely. For example, corporate 
asset correlations, which determine the correlations of stock and corporate bond returns, 
are usually positive. This is partly for economic reasons but it also has an underlying 
mathematical origin. If a N x N correlation matrix has all its off-diagonal elements equal 
to £ then it is only mathematically consistent (i.e, has non-negative eigenvalues) for — l/(iV — 
1) < £ < 1 [Wise and Bhansali (2002)]. In the limit N — » oo negative correlations are not 
allowed and hence the average correlation (found by averaging over allowed values of £ with 
any nonsingular probability distribution) will not be zero. 
The expected value of P's relative return vanishes, 



E[AR] 



0. 
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However Vs relative return does fluctuate. Define the moments of P's relative return 
probability distribution by, 
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The second moment, m = 2 is the square of the tracking error. For m = 2 [Dynkin, Hyman 
and Konstantinovsky (2002)] , 
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2 r (2) - g re i a ted to the covariance matrix p by, pij = 
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For the coefficient of r\ 2 the first term in the square brackets comes from terms in the sum 
of equation ([?]) with both the subscripts of corresponding to assets in the portfolio V, 
the second term comes from terms in the sum of equation (^) with one of the subscripts 
of r^ 2 - 1 corresponding to assets in the portfolio V and the third term comes from terms in 
the sum of equation ([?]) with none of the subscripts of corresponding to assets in the 
portfolio V. 

Taking the limit TV — > oo gives, 

«< 3 > = ±(rg> - rg>). (9) 

Even with non-zero asset return correlations diversification achieved by increasing n reduces 
the variance of Vs return relative to the benchmark. This is very different from the behavior 
of the variance of Vs return (without the benchmark subtracted) which, in the presence of 
significant correlations, cannot be made small by increasing n. 

We do not assume that the probability distribution for Ai2 is normal and so the moments 
t;( m ) with m > 2 are not related to the tracking error. The purpose of this paper is examine 
the behavior of these moments of the probability distribution for Vs relative return when 
N » 1, n >> 1 and n/N « 1. This is particularly important in times of market stress 
for portfolios containing corporate bonds. In such eras default correlations are large [Das, 
Freed, Geng and Kapadia (2001)] and as remarked earlier the joint probability distribution 
of corporate bond returns is expected to be very far from multivariate normal. 

We are interested in the value of the m > 2 moments of Vs relative return as the 
benchmark portfolio gets arbitrarily large, N —* oo, and write in that limit 

(m) i \ (m) i / \ (m) . , \ (m) 

+ a 2 ,2(m,n)r[™ 2 ) 23 _ m _ 2 + . . . . (10) 

The subscript on a denotes how many subscripts of take on the same value. For 
example c^...^ corresponds to i\ indices taking on the same value (which without loss of 
generality can be taken to be 1), ii of the indices taking on the same value (which without 
loss of generality can be taken to be 2), etc. Equation @ yields, ao(2, n) = — 1/n and 
0:2(2, n) = 1/n. Performing a similar computation for the third moment yields, ocq(3, n) = 
2/n 2 , 0:2(3, n) = —3/n 2 and 0:3(3, n) = 1/n 2 . 

Generalizing the N —> 00 limit of the coefficient of r^ 2 in equation (|8|) to arbitrary m 
gives (assuming n > m), 

« = g(^^)(^) (-!)-'(! " 

= (-l) m n- m U(-m,l-m + n,n), (11) 

where U is the second solution to the confluent hyper geometric equation^]. The k = p term 
in the sum corresponds to terms in the sum of equation (0) where p of the subscripts of 
r h^i m l a b e l assets in the portfolio V. The remaining m — p subscripts label assets not in 



See, for example, pages 753-758 of Arfken (1985). 
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V . With m held fixed and n taken large U(—m, 1 — m + n, n) is of order rS m ~ 1 ^ 2 for m 
odd and of order n m / 2 for m even. 
This suggests that for n » 1, 

(2m) „(2m+l) j 

(» ,2> ) ~ OWn). ^ ~ 0(1), and ~ O^). (12) 



and it is not difficult to show this is indeed true by also considering the cases where some 
of the subscripts on r| 4 take on the same value. If j indices take on the same value then 
to get a contribution that survives in the limit N — > oo these indices must correspond to 
assets in V and so in this case, 

/ m \ \ ^p 7 / ( m — j)\ \ / (n — 1)1 \ _■ fc/lV _1+fc 

aj(m ' n) = £Ul(m-j-fc)!j l(n - 1 - fc)! J ' + U 
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(-l) m -V~ m £/(j -m,j-m + n, n). (13) 



Since for large n the function f7(j — m,j — m + n, n) is of order n*-" 1 J for m — j odd 
and of order n( m- j)/ 2 for m — j even these terms leave unchanged the behavior in equation 
dH). The other cases are handled similarly. 

A simple example illustrates the power of the constraint on the probability distribution 
for Vs relative return that equation (O) provides. Consider the probability distribution, 



P(AR) 
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where 5 denotes the Dirac delta function and a is a positive constant. For large n the second 
moment of this probability distribution is ~ (1 + a) /(An) and for m > 1 (and n large) 
its (2m)'th moment is r/ 2m ) ~ a/(2 2m n). While this probability distribution is consistent 
with the behavior the second moment in equation ( |l^ ) for a ^ it is in conflict with the 
behavior of the higher moments. With a / the probability distribution in equation ( |l4| ) 
implies that the probability for the relative return AR = ±1/2 is of order 1/n, which is too 
large to be consistent with the scaling of the higher moments in equation (|T^). 

We have generalized the results of Dynkin, Hyman and Konstantinovsky (2002), estab- 
lishing that higher moments of Vs relative return become successively smaller for n >> 1 
(by powers of 1/i/n) even if asset return correlations are large and the joint probability 
distribution for the asset returns is not multivariate normal. From equation (|T^) it is also 
evident that as n gets large the skewness becomes small (i.e., 0(1/ y/n)). For N » 1, 
n >> 1 and n/N « 1, the probability distribution for Vs relative return is approximately 
symmetric and strongly peaked about zero. However, it does not necessarily approach a 
normal one (e.g. the excess kurtosis need not be small). For very large n the scaling 
in equation ( |l2| ) is consistent with the probability distribution for AR being a symmetric 
function of y/nAR. 



4 



References 



Arfken, G. (1985) Mathematical Methods For Physicists, Academic Press Inc. 

Black, F., and Cox, J. (1976) Valuing Corporate Securities: Some Effects of Bond Indenture 
Provisions, Journal of Finance, 31, 351-367. 

Das, S., Freed, L., Geng, G. and Kapadia, N. (2001) Correlated Default Risk, working paper, 
Department of Finance Santa Clara University and Gifford Fong Associates. 

Dynkin, L., Hyman, J. and Konstantinovsky, V. (2002) Sufficient Diversification in Credit 
Portfolios, Lehman Brothers Fixed Income Research. 

Leland, H. and Toft, K. (1996) Optimal Capital Structure, Endogenous Bankruptcy and the 
Term Structure of Credit Spreads, Journal of Finance, 51, 987-1019. 

Longstaff, F. and Schwartz, E. (1995) A Simple Approach to Valuing Risky Floating Rate 
Debt, Journal of Finance, 50, 789-819. 

Wise, M. and Bhansali, V. (2002) Portfolio Allocation To Corporate Bonds with Correlated 
Defaults, To appear in the Journal of Risk. 



5 



